Transport of exotic anti-nuclei: II- p and d astrophysical uncertainties 
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^ Abstract 

^ I We use a ID propagation model to study the dependence of the p and d exotic fluxes on the transport parameters. The simple 
— analytical solutions allow us i) to clarify the origin of the astrophysical uncertainties, and ii) to compare two models used for 
fSJ signal predictions, namely the constant and the linear Galactic wind models. We also study how these uncertainties should be 
T-H reduced using forthcoming nuclear cosmic ray data. We confirm that the degeneracy of the transport parameters for a given 
propagation model leads to very different fluxes for primary antinuclei (--^ 10^). However, we show that with forthcoming data, 
these uncertainties could be greatly reduced (~ 2). As the precision will increase, the astrophysical uncertainty could then be 
(-H nominated by our ignorance of the correct spatial dependence for some of the transport parameters: for instance, the constant and 
^ ^the linear wind models do not predict the same amount of exotic p at low energy. 
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1. Introduction 

' Positron, anti-proton and anti-deuteron fluxes are among 
the primary targets of ongoing (bess, Pamela) and forth- 
coming experiments (ams, gaps), seeking for indirect liints 
at dark matter candidates. These particles, unlike 7-rays, 
are very sensitive to the random magnetic fields pervading 
the Galaxy, and their transport is diffusive in nature. 

Several propagation models (e.g., DarkSusy — Gondolo et 
2004; Galprop — Moskalcnko et al. 2005; or Donato ct al. 
200f ) are being used to calculate the standard background 
fluxes as well as possible exotic contributions from new 
particles filling the dark matter halo. These calculations 
aim to uncover excesses over the background, or to put 
some constraints on parameters of new physics if none 
is observed. In this context, it is crucial to understand 
the uncertainties associated with the background and the 
hypothetical signal. Throughout this paper, we focus on 
astrophysical uncertainties only. To estimate them, we 
assume that the source spectrum and spatial distribution 



of the exotic matter are fixed, leaving the propagation 
coefficients as the only free parameters of flux calculations. 

Whatever the propagation model, the exotic fluxes (and 
astrophysical uncertainties) depend on the associated 
transport coefficients (and their uncertainties). The usual 
way to derive all these quantities is as follows: 

1. Transport parameters are determined using secondary 
to primary ratios, usually B/C, i.e. from standard GCRs 

g^l whose sources are located in the disk of the Galaxy. A 
numerical study is generally required to extract 

(a) the parameters yielding the best fit to data; 

(b) uncertainties on these parameters. 

2. These parameters are then used to calculate 

(a) e+ , p, and d secondary fluxes (background) , whose 
sources are also in the Galactic disk; 

(b) e'^ , p, and d exotic fluxes ^ , whose sources are now 
located in the diffusive halo of the Galaxy. 

Transport coefficients from the best fit give the most 
likely (or p d) fluxes, whereas the range of allowed 
transport parameters determine their uncertainties (for 
the chosen modelling) . 
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Below, we use indifferently 'exotic' or 'primary' fluxes. 



A thorough study along the hues sketched above has 
been performed in the framework of the 2D constant wind 
diffusion model (Maurin ct al., 2001; Donato ct al., 2001; 
Barrau et al., 2005). With respect to item (1), Maurin ct al. 
(2001) found that the set of allowed propagation parame- 
ters is strongly degenerated. As regards (2a), Donato ct al. 
(2001) found that this degeneracy has almost no conse- 
quence on the secondary p (or d) fluxes (uncertainties are of 
about 25% at most) — which is expected for any modelling, 
as long as the propagation parameters are adjusted from 
B/C data. Finally, related to (2b), Barrau ct al. (2005) 
showed that this degeneracy leads to an uncertainty on ex- 
otic p fluxes of about two orders of magnitude at low energy. 
The purpose of this paper is partly to explicit with simple 
calculations the above result, then to extend the analysis 
to different propagation models. 

The outline of the paper is the following: in Sec. 2, we 
start with a description of the ID modelling used through- 
out the paper. Section 3 deals with astrophysical uncer- 
tainties of p and 0? primary exotic fluxes. The analysis is 
performed for the case of pure diffusion (§3.1), then for the 
case of diffusion/convection models (§3.2). The predictions 
for the constant and the linear wind models are compared 
in Sec. 4, using the best-fit values of the transport param- 
eters provided in the literature. Finally, in Sec. 5, we focus 
on the constant wind model and estimate by how much the 
existing uncertainties on the p and d exotic fluxes could be 
reduced with better data. 

2. ID diffusion models 

Diffusion models in which the density is assumed to de- 
pend only on the z-coordinatc (ID models) provide a sim- 
ple but useful description of the propagation of cosmic ray 
nuclei (Jones, 1979; Jones ct al., 2001). 

Regarding exotic sources, we showed in the companion 
paper (Maiu'in ct al., 200G) that analytical ID-models pro- 
vide a fair description of the p primary fluxes, compared to 
that calculated from a full 2D-modcl taking into account 
energy gain and losses. The main reason is that the energy 
redistribution terms arc restricted to the thin disk. How- 
ever, exotic sources of p arc located in the extended diffu- 
sive halo of the Galaxy, and the p reaching us rarely cross 
the disk on average. 

Thus, if we discard energy gains and losses and use a ID 
geometry, the p and d flux calculations will still be satis- 
factory for a gross estimate (compared with full 2D calcu- 
lation), but now the advantage is that we deal with simple 
analytical solutions (sec App. B and C). These solutions al- 
low us to put to the fore the key parameters driving exotic 
anti-nuclei uncertainties, as well as a their quick estimates. 

2.1. Short description 

We consider an infinite plane whose density and source 
distribution do not depend on r (see Fig. 1). The thin-disk 
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Fig. 1. Upper half-plane of the ID infinite disk along r (half-thickness 
h). Cosmic rays diffuse in the disk and in the halo, with a difliusion 
coefficient K(z,E). 

approximation is used (h <ti L). The diffusion/convection 
transport equation (without energy redistributions) for a 
constant source term goark (E) reads 

-^[Ki^)^]+^[Vs.iiz)N]+2hrd{z)N = q^,,^{E)il) 

The spatial and energy dependence of the diffusion coeffi- 
cient are assumed to be independent: K{z, E) ~ I3Kq{z)TZ^ 
{TZ is the rigidity, f3 = v/c). 

Discussions in the paper will refer to the four propagation 
parameters of the model: i) the halo size of the Galaxy L, ii) 
the diffusion coefficient normalization Ko(z), iii) the diffu- 
sion slope 6, and iv) the value of the Galactic wind Vgai{z). 
Many propagation models arc equivalent to a Leaky Box 
Model (LBM), as long as stable species originating from the 
disk are considered (Jones, 1970). Hence, transport param- 
eters of ID-models can be easily inferred (step 1 described 
in Sec. 1) from the mean grammage {x)^^ = Xcsc{E) of 
the LBM: the details of this equivalence as well as specific 
relationships are reminded in App. A; mean grammage in 
various ID configurations are given in App. B. Step (2) of 
the analysis, i.e. estimating uncertainties on exotic fluxes, is 
then straightforward using ID formulae derived in App. C. 

For the sake of simplicity, spallative destruction of exotic 
species in the disk is neglected throughout the paper. It 
is checked in App. D, for the constant wind model, that 
this leaves the conclusions unchanged. This is understood 
as the bulk of exotic p reaching us are mostly those created 
in the halo: on average, their disk crossings are too rare to 
yield a sizeable effect on the low energy spectrum (energy 
redistribution are neglected for the same reason) . 

2.2. Type I and Type II uncertainties 

Two types of uncertainties are identified and discussed 
in this paper: 

Type I (parameter uncertainty) corresponds to uncertain- 
ties related to the degeneracy of transport parameters for 
a given propagation model (i.e. departure from best fit). 
Type II (modelling uncertainty) corresponds to the differ- 
ent predictions of primary fluxes between different prop- 
agation models whose respective transport parameters 
give the same standard ffuxes. 
We focus on the following conflgurations: pure diffusion 
(with one diffusion coefficient) and constant or linear wind 
model. 



3. Degeneracy in the transport parameters 

The diffusion coefficient is assumed independent of the 
spatial coordinates, so tliat Kq{z) = Kq. 



One of tlie main appeal for adding a constant wind comes 
from the measured form of Aosc(^')- A fit to the B/C data in 
the Leaky Box Model — including all ingredients, i.e. energy 
losses, etc. — gives (see, e.g., Webber et al. 1998) 



3.1. Pure diffusion 

3.1.1. Degeneracy of parameters for standard sources 
In the thin disk model, the transport equation reads: 

~KN" + nva2h5{z) x N = 2hS{z)q. 

For z = 0, the solution is 



iV(0) = 



2hq 



2K/L + nva ' 

Inserted in Eq. (A. 6), this gives for the mean grammage 

/ \ pure- DM _ f^^'^ 



2K 



(2) 



where ^ = 2hnm is the surface mass density in the disk, L 
the halo size, K the diffusion coefficient and v the velocity 
of the nucleus. 

As the mean grammage in the LB model is given by 
(x)^^ = Acsc(£^), the diffusion coefficient is inferred from 
Eq. (A.7)— i.e. {x)p 



ure-DM ^ l^x)^B (2): 
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(3) 



2\csc{E) 

Both the diffusion model and the leaky box model lead to 
the same B/C ratio as long as Eq. (3) holds. Consequently, 
both models lead to the same standard secondary fluxes for 
p and d. 

3.1.2. Lifting the degeneracy: uncertainties for primary p 
and d in the halo 

In contrast, exotic sources are located in the whole diffu- 
sive halo, which lifts the above mentioned degeneracy. The 
equation (no spallations) now reads 

-KN" = goark, 

whose solution is given by Eq. (C.8) 
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For a given value of Acsc, i-c a given B/C ratio, the pri- 
mary ffux depends linearly ^ on L. If a very conservative 
range, 1 — 15 kpc, is assumed for L (as taken for example 
in Maurin ct al. 2001), then the primary flux is uncertain 
by a factor of 15. 

3.2. Adding a constant wind Vga\ = Vc 

Keeping the same spatial dependence for the diffusion 
coefficient, a constant wind directed outward the Galaxy is 
now added (pure diffusion recovered when Vc — > 0) . 



This behaviour can be seen in Fig. 8 of Bai i au ct al. (2002) for a 2D 
model using a realistic dark matter distribution (for PBH sources). 



/3 X 15.61 g cm-2 if7^<3.6GV; 
Aesc(7^) = { (5) 



P{TZ/3.6) "'^x 15.61 g cm '^otherwise. 

Such a form naturally arises in diffusion/convection mod- 
els (compare (x)^" and Xcsc{E) at low energy). This re- 
markable result was first pointed out by .Jones (1979): the 
phenomenological quantity Acsc(^), taken from the LBM, 
hints at a dual transport in propagation models (see also 
Jones et al. 2001 for other interpretations). Beyond a few 
GV, the shape of Aesc shows that diffusion always prevails 
so that conclusions of Sec. 3.1 still hold at high energy. 
The mean grammage is obtained from Eq. (B.8) 



^ ' - 2Vr 



1 - e" 



(6) 



This formula has two asymptotic behaviours with energy 
{K{E) depends on energy): i) a diffusion-dominated regime 
\iVcL/K <C 1 (reached at high energy) and ii) a convection- 
dominated regime if V^L/JC 3> 1 (low energy) . The limiting 
cases are 



fivL 
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and 



(7) 



At high energy, the same interplay between the propa- 
gation parameters applies as in Sec. 3.1. At low energy, 
an upper limit for the wind velocity is found when the 
convection-dominated regime is reached: equating rhs of 
Eq. (7) and first line of Eq. (5) , with a standard value ^ = 
g cm~^ gives. 



2.4 X 10-3 " 



Vc - 20 km s" 



(8) 



Equation (C.8) gives the fiux for exotic sources in the 
galactic halo (as before, spallations neglected): 



AT^^ (z^O) '^D"'^^'' 
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which asymptotic forms read 
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(10) 



Equation (6) gives VcL/K as a function of Aosc, provided 
that V < IJ'V /2\csc (^ 20kms-^). Equation (9) is rewritten 
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(11) 



For a fixed value of Acsc, the exotic fiux depends linearly on 
K and only weakly on the combination VcL/K. This gives 
an overall uncertainty similar to the asymptotic case. 



3.3. Linear wind Vj,ai(^) = VJ x z 

A linear wind Vgai(z) = V; x z is now considered 
([V;] =km s~^ kpc~^). The mean grammage is obtained 
from Eq. (B.9) 



^jLVy/iT erf 



V 2K 



and the primary flux is given by Eq. (C.7) 



The corresponding hmiting cases are 



2K 



v,«- 



'2K 



and 



v,»- 



V8KV1 ' 



(12) 



(13) 



(14) 



and 
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At low energy, we proceed as in the previous section (at 
high energy, the pure diffusive transport is recovered). It 
is first assumed that the conveetive transport regime Vi ^ 
K/LP' holds and that the mean-grammage is perfectly de- 
termined from hypothetical perfect data. The mean gram- 
mage fixes Vi through right-hand side Eq. (14). This time, 
the degeneracy on Kt^/ L leads to no uncertainty on the ex- 
otic flux, at variance with the former case: 
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compared to 

Vi 
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qK 



This simply pertains to the absence of the factor L or if in 
the left-hand side formula. 

As in Sec. 5.1, we can derive estimates of the range span 
by Vi from real data. At low energy, one should obtain from 
Eqs. (5) and (14) the relation 

15.61 g cm-2 ^ 4^ . 

From the range of variation of the LB model (i.e. 5 = 0.3 — 
0.7) using /i = 2.4 x 10~^ g cm~^, we get 

Vi -8-30 kms"^ kpc-^ 

This value is not too sensitive to the uncertainty in L be- 
cause the latter translates in Vi through \fK in the mean- 
grammage formula above. This is compatible with results 
gathered in Tab. 2 of Lionetto et al. (2005) (for 8 0.3). 
The "parameter" uncertainties (Type I) for this model are 
moderate (a factor of 3 here) and smaller than at high en- 
ergy: it was the reverse for the constant wind model. 

4. Difference between models: constant vs linear 
wind 

The constant wind and the linear wind models are 
both used to exclude (or fit) parameters of new physics 



(Donato et al., 2004; Lionetto et al., 2005). It is thus of 
paramount importance to check if they predict the same 
amount of exotic p and d. 

This difference exists only at low energy when the two 
models differ from pure diffusion. Taking the ratio of the 
primary flux formula of the linear-wind model Eq. (13) to 
that of the constant-wind model Eq. (9) leads to 
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where Ki and Kc are the value of the diffusion coefficients 
found respectively in the linear and constant wind model. 
As a new exotic component is usually seeked at low energy, 
we only calculate this quantity for a kinetic energy Ek ~ 
600 MeV (IS), corresponding roughly to 7?. — 1 GV, so that 
K{1 GV) - Kq. 

The diffusion coefficient K should be constrained to be 
the same in both models (high energy constraint). Table 2 
of Lionetto et al. (2005) gives A'o = 2.5 x 10^^ cm^ s"\ 
L = 4 kpc and Vi = 6 km s~^ kpc~^ (for 6 = 0.55) whereas 
for the constant wind model. Tab. II of Barrau et al. (2005) 
gives A'o = 1.35 x 10^^ cm^ s"\ L = 4 kpc and Vc = 
12 km s^^ (for 5 = 0.7). Note that both best fits correspond 
to L = 4 kpc (although it could be only a coincidence). 
This minimizes the difference that is observed between the 
two calculations as different L give different primary fluxes 
(see Sec. 3. 1.2). Plugging these parameters in the previous 
formula gives 



(P, d] 



^\p, d) 



1.3 . 



This is not striking a difference, but this figure is expected 
to increase at lower energy. Moreover, it would also be larger 
would the best-fit value for L in the two models be different. 

In any case, we remind that in any models, the transport 
parameters are more or less degenerate, providing the right 
amount of B/C. However, as explained in a previous sec- 
tion, this degeneracy is broken for (exotic) sources in the 
Galactic halo. The way it is broken depends on the vertical 
dependence of the transport parameters. Hence, different 
models are expected to lead to different predictions for the 
exotic flux. This issue deserves a dedicated study, which 
goes beyond the scope of this first simple discussion. 

5. Reducing uncertainties 

5.1. Impact of the inaccuracy of B/C data 

So far, we discussed the hypothetical case of ideal data 
to estimate the degeneracy in the transport parameters. 
Using real ones, a greater degeneracy is generally obtained, 
enhancing the uncertainties on the exotic fiuxes. 

In LB models, the acceptable range of parameters found 
in the literature providing a good fit to B/C is i5 = 0.3 — 0.7. 
This value of the spectral index in Aesc(''^) [see Eq. (5)] 
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0.85 
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Table 1 

Propagation parameters consistent with B/C data (Maurin ct al., 
2001). The set labelled best corresponds to the best fit to B/C data, 
while those labelled min and max correspond to sets which give 
minimum and maximum exotic fluxes (Donato ct al., 2004). 



read from Fig. 2, reaching a maximum at low energy and 
decreasing at high energy. For comparison purpose, uncer- 
tainties from the 2D modelhng are also displayed (using the 
same propagation parameters; see details about the 2D- 
model in the companion paper) . For standard dark matter 
profiles, the two approaches are in fair agreement . The 
various numbers found strengthen the use of the ID-model 
for qualitative estimates of uncertainties. 

5.2. How to decrease the astrophysical uncertainties? 
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Fig. 2. Ratio R of the primary fluxes calculated from the max and 
the min propagation parameter sets. The solid line results from 
the ID modelling Eq. (9), while all other curves come from full 
2D calculations using three different dark matter profiles (adapted 
from (Barrau ct al., 2005)). 

depends on the level of reacceleration (.Jones ct al., 2001). 
This translates into an increase of the range spanned by 
Aosc (??•), and thus to an increase of the propagation pa- 
rameter range extracted from it in the ID model. We do 
not pursue along this line. Rather, taking advantage of the 
1D«-^2D equivalence, we directly read these propagation 
parameters from the 2D-model, as given in Maiuin ct al. 
(2001). For reasons exposed in Sec. 3.2, configurations of 
maximum wind strength Vc and minimum halo size L lead 
to the minimal exotic flux and vice versa. Table 1 gathers 
the set of parameters leading to the maximum/minimum 
exotic fluxes, as well as the configuration corresponding to 
the best fit to B/C. 

Using these max and min configurations in Eq. (9) 
for the low energy primary exotic flux leads to the result 
N {max) / N [min) = 90. Compared to w 10 which was ob- 
tained from the pure-diffusive regime, this implies that Vc 
increases the uncertainty on primary fluxes when S is not 
well constrained . 

This ratio can be estimated at all energies. The corre- 
sponding "energy-dependent" uncertainties can be directly 



^ We underline that this is in fair agreement with the results ob- 
tained i) in a semi-analytical similar — albeit more complicated — 
study within a 2D propagation model or ii) using the full propaga- 
tion scheme in this same 2D model (see respectively Sec. III.C and 
Sec. IV in (Donato et al., 2004)). 



The above result is a further support to our approach. 
It is thus tempting to pursue along this line and estimate 
how the uncertainties could be decreased with better data. 

5.2.1. Sensitivity to the diffusion slope 5 

The uncertainty on the value of S greatly enlarges uncer- 
tainties on exotic fluxes at low energy. Castellina and Donato 
(2005) showed that the ongoing CREAM experiment (Seo et al. 
2004) should pinpoint 6 with a precision of 10 — 15%. 
Reading the corresponding intervals span by the transport 
parameters from Figs. 7 and 8 of Maurin et al. (2001), a 
simple calculation — analog to what is done in the previous 
subsection — shows that the uncertainty at low energy is 
decreased from the previous value of ~ 100 down to ~ 10. 
The latter value is in agreement with the estimate for 
fixed S, meaning that the only remaining uncertainty 
comes now from the Kq/L degeneracy. The PAMELA ex- 
periment (Picozza ct al., 2000), which was successfully 
launched in June 2006, or AMS on ISSA could as well help 
determining with a very good precision this parameter S. 

5.2.2. Radioactive nuclei 

It is well known that radioactive nuclei, e.g. the ^"Be/^Be 
ratio lifts the K/L degeneracy. To complete the exercise, let 
us investigate to which precision L could be obtained. To 
date, the best data come from ACE (Yanasak et al., 2001): 
they correspond to a few energy points only — as it is the 
case for almost all radioactive data. We use the ^"Be/^Be 
ratio, as heavier isotopes (such as the ^^Al/^^Al ratio) suffer 
larger experimental uncertainties. 

Some quick procedure was exposed in .Jones (1979) to 
extract the value of Kq (or L) and its uncertainties from 
the data. In our case, as we are only interested in the ratio 
of extreme values, i.e. L'^'^^/L™™, the discussion is even 
simpler. What is demanded is that 



10 



Be 



'Be 



10 



Be 



'Be 



ID-ii 



3dcl 



Expressions for the radioactive species '^Be and the stable 
isotope ^Be remain to be determined in the ID framework. 
For ^°Be, as the decay time jtq is generally much smaller 



* The minimum in the curves appear because propagation models 
were demanded to fit the 10 GeV/n B/C data point. 



than the convcctivc, spallativc or escape time, the equation 
to solve in the ID-model is simply 



K- 



+ jtqN = 2hqsec5{z) 



whose solution is 

^rad.n^ 



hq 



, tanhfL.O^V- , . 



(16) 



where the approximation comes from L^^tq/ K > 10, so 
that tanh(. . . ) « 1. 

For the associated stable isotope, the solution can be 
found in the Appendices. However, it is not even needed: 
^Be is a stable secondary nuclei having the same propaga- 
tion history as boron. Demanding that a model fits B/C — 
which is a pre-requisite at this stage — automatically en- 
sures the constancy of the ^Be flux. It means that 
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Expressing K = KqPR^ and getting the constant ratio 
KqIL (fitting B/C) to appear in Eq. (16), we obtain 
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Using now for the ^"Be data the 3 — cr error bars given 
in Yanasak et al. (2001) results in 



0.120 + 3 X 0.008 
0.120 - 3 X 0.008 



2.25 



Hence an uncertainty of about a factor two in L (or equiva- 
lently in Kq), and the same uncertainty in the exotic spec- 
tra. 

However, it should be kept in mind that such a small 
uncertainty is correct as long as b is known. For instance, 
Donato ct al. (2002), using their derived range of allowed 
values for (5, found that the use of radioactive data could 
not lift the A'o jL degeneracy (although it slightly reduces 
it). Moreover, there could be other issues with radioactive 
data, as for example the effect of specific /ocaZ features in the 
solar neighbourhood (see, e.g., Donato et al. 2002), which 
could affect the determination of the transport parameters. 

This exercise has been conducted using existing ACE 
data. Future instruments, such as AMS have the capability 
to achieve better performances. 

5.2.3. Summary for the constant wind model 

Table 2 summarizes the qualitative results derived in this 
section. Determining the slope of the diffusion coefficient 
from B/C data should be the primary concern, whereas ra- 
dioactive nuclei ultimately lift the remaining degeneracies. 

We also remind that, although it was not explicitly stated 
all along, all the reasonings do hold for exotic anti-protons 
as well as for exotic anti-deuterons. 
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changed whether radioactive nuclei arc considered or not. 

^ Astrophysical uncertainty should vanish using new radioactive data. 

Table 2 

Summary of results derived in Sees. 3.1 and 3.2. The first row cor- 
responds to (1) left column: formula for the grammage fitted from 
B/C data and (2) right column: the primary fiux when spallative 
destruction is discarded. After the second row, each column (1 and 
2) is further subdivided in two parts corresponding to the asymp- 
totic behavior i) at low energy (left; convection-dominated) and ii) at 
high energy (right; diffusion-dominated). Note that the convection- 
dominated regime is not necessarily reached. In the subsequent rows, 
col. (1) corresponds to the degeneracies (derived from the above for- 
mulae) in the transport parameters and col. (2) to the uncertainties 
associated with these degeneracies. The same structure (further sep- 
aration into low and high energy) is kept as for the upper half of 
the table. For example, for current B/C data at high energy, the de- 
generacy is fi ■ L/Ko, leading to a relative uncertainty of AL/L for 
exotic species. 

6. Summary and conclusions 

We used ID-propagation models to study astrophysical 
uncertainties on exotic p and d fluxes (e.g. from annihilation 
of dark matter in the Galactic halo). We started with a 
detailed analysis of the constant wind diffusion model to 
recall and understand the origin and key parameters of the 
uncertainties for exotic fluxes. 

We then stressed the origin of uncertainties for the 
transport parameters in the constant and linear wind 
models, and their differences at low energy, emphasizing 
on how crucial the determination of the diffusion slope 
S was. We showed that the present astrophysical uncer- 
tainties of 0(100) could be severely reduced with a good 
determination of S. 

We finally showed that even if the internal degeneracies of 



a given model are reduced, the differences between different 
modellings (e.g. assuming different spatial dependences of 
some of the transport coefficients) could be important. For 
the constant and the linear wind model, a lower limit of 30% 
was derived. These differences between models are bound 
to become a more serious matter in the near future. 



measurements are made, i.e. Tq. Thus, a diffusion models 
(DM) will be locally equivalent to a LBM if 

G^^^(x,r0) = G^^*^(a;) . (A.4) 

We remind that for a LBM all locations come to the same 
thing since all quantities are spatially averaged. 



A. 2. PLD and mean grammage (x) 
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Appendix A. Reminder of the equivalence 
between Leaky Box and diffusion models 

The link between diffusion models and Leaky Box Mod- 
els is established through the weighted-slab formalism. The 
weighted-slab decomposition is only approximate when en- 
ergy losses are taken into account (Lczniak, 1979; .Jones, 
1991). But, in this paper, this is not an issue as they are 
neglected for all ID models. 

A.l. The weighted-slab formalism in brief 

The ansatz is to rewrite the differential density N{r) 
of any transport equation (e.g. Eq. 1) as N{r) = 
Jp N{x)G{x,r)dx (where x is the grammage crossed by 
a CR in g cm~^). Demanding N{x) to follow the slab 
equation, one is left with an equation on G{x,r). The 
slab equation is common to all models (LB, ID and 2D- 
diffusion models, etc.) and describes the nuclear aspects 
of propagation. The remaining equation only depends on 
the transport properties and on the choice of the geome- 
try, both assumed to be independent of the species (see 
Sec. 2 of Mam-in ct al. 2001, for a longer albeit succinct 
summary; or Berezinskii et al. 1990, p. 45 for a more com- 
plete treatment). The exercise is easily done on the leaky 
box equation (no spatial dependence) governing the evo- 
lution of the nucleus j (source term plus secondary 
contributions from all heavier nuclei k) 



^-^ m 

k>] 



(A.l) 



Using the ansatz and inserting the slab equations 



dN^jx) 
dx 



kj 



-N^x)^qi+ y —N''{x) (A.2) 



k—j-\-i 



leads to an equation for G{x) whose solution is given by 



G^^*'(^) = T^cxp(-^ 



(A.3) 



Different propagation models can be compared through 
their path length distributions at the location where CR 



A useful quantity is the mean grammage crossed by CRs, 
defined as 

xG{x, r)dx 



G[x,r)dx 

An alternative expression is provided by 
d 



(x)(r) 



daP 



IniVP(r) 



(A.5) 



(A.6) 



aP=0 



where NP{r) is the solution for a primary species (see 
Berezinskii et al. 1990, p. 53). The advantage of this sec- 
ond formulation is that there is no need to know the PLD, 
G{x,r), to evaluate the mean grammage. Indeed, the for- 
mer quantity can be trickier than NP{r) to derive (as an 
illustration, PLDs for several modellings can be found, e.g., 
in Owens (1976)). 

The mean grammage is only the first moment of the dis- 
tribution G{x,r): only the latter contains the full infor- 
mation on the transport for the studied model. However, 
for two propagation models to be equivalent, we will use 
a weaker version of Eq. (A.4), only demanding the equiva- 
lence of the first moments 



\DM 



\LBM _ 



= Acsc(-E) 



(A.7) 



The last equality comes from Eq. (A.3) plugged into 
Eq. (A.5). We remind that all the above quantities de- 
pend on the energy. Note also that this relation does 
imply equivalence of PLDs for many cases especially if the 
thin-disk approximation is assumed (see App. B.2). 

Appendix B. Mean grammage (x) in two 
modellings 

The diffusion coefficient K is the same in the disk and in 
the halo (see Fig. 1) . The galactic wind is taken either as 



^ Contrarily to an erroneous calculation in tiie previous versions 
of this paper, in the thin disk limit, setting two distinct diffusion 
coefficients for the disk and for the halo is strictly equivalent to have 
only one diffusion coefScient for the two zones. Indeed, in the thin 
disk approximation, the time spent in the disk tends to zero, while 
the number of crossing (hence the grammage) is kept constant: only 
the diffusion coefficient in the halo plays a role. This result is of 
course known for almost as long as diffusion models exist. The trivial 
mistake we made (shame on us... and oh, yes, a referee to point out 
the mistake would have been just the right thing to have, but as 
refering goes, you know...) was connected to a wrong implementation 
of the continuity condition at the crossing of the zones. A very small 
part of the paper was affected by this mistake, and the only one 
plain wrong conclusion is now removed. 
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a constant Vgai^z) ~ Vc or as a linear term Igai(^) = V; x z. 
The transport equation reads: 

-KN" + (VgaiiV)' + 2h5{z).nv(j x N = 2hqS(z). (B.l) 

The quantity nav is the total destruction rate of the anti- 
nucleus under study with the interstellar gas, whose density 
is n. 

B.l. Solutions 

For our purpose, it is sufficient to extract solutions for 
primary standard species, that lead straightforwardly to 
the mean grammage (a;) . 

B.1.1. Constant wind Vgia{z) = Vc 

We start with the equation in the galactic halo, i.e. 

-KN" + VcN' = N = Ae^'"-^^^ + B . 

Implementing the boundary condition iV(z = i) = de- 
termines one of the two constants of the solution in the 
halo: 

N^A f e^=^/^ - e^^^/^ 



We integrate over the thin-disk (see e.g. (Jones et al., 
2001)), i.e. lime^o /J^/(Eq. ^.\)dz, taking care of the dis- 
continuity induced by the wind during the disk crossing: 

-~2KN\^) -t- 2VcA^(0) H- 2hnvaN{}S) = 2hq . 

The solution for the halo N above is plugged into this equa- 
tion. It is convenient for the following to write the solution 
as: 

l-p--^(i-^) 

N{z) = A^(G) X 



1 - e- 



with 



NiO) 



2hq 



2Vc/(l - e-^-i/-^) + 2hnva 



(B.2) 



(B.3) 



B.l. 2. Linear wind Vgai(2:) = Vi x z 

The parameter Vi has the dimension of the inverse of a 
time. As above, we start by solving this equation in the 
halo, where it reduces to 

-KN" + VizN' + ViN = . 

This equation is not hermitian and it is more convenient to 
use new variables n and y defined as 

n{y) — N{y) and y = \/2vz 

where 

The diffusion equation then reads 



(B.4) 



n"(?/) - n{y) 



4 



. 



The constants of integration A and B are determined as 
usual. First, the condition N{z = L) = Q yields^ 

A = -B X erf (i^L). 

Second, integration of Eq. (B.l) over the thin disk yields 

-2A'iV'(0) + 2hrN{0) = 2hq , 

which determines A and gives the final expression 



N{z) = N{0) X e^"^^^ 



with 



iV(0) 



eri{vL) — erf(i^z) 
erf (jyX) 



2hq 



where f is defined in Eq. (B.4). 



-I- 2hnva 



(B.5) 



(B.6) 



B.2. Mean- grammage [x] at z = 

Two approaches can be used to derive the mean gram- 
mage. The first one is depicted in App. A. The second one 
is based on the fact that, in the thin disk approximation, 
all solutions for z = can be recast to look formally as to 
a Leaky Box solution (see .lories ct al. 2001): 

2hq/ 



N{0) 



ID-r 



odd 



(B.7) 



l/(a;)iD-modoi _^ ■ 

In that case, demanding an equivalence between the 
mean grammages of different models (as in App. A. 2) is one 
and the same to demanding the equivalence between path 
length distributions. The corresponding mean-grammages 
are directly obtained from Eqs. (B.3) and (B.6). Using the 
relation fi = 2hnfh (fh is the mean mass of the interstellar 
medium and fj, is the surface mass density of the gas in the 
disk), 

- Constant galactic wind Vga\{z) — Vc 



^ ^ 214 V 
Linear galactic wind Vgia{z) = ViXz 

fiv X V^erf (^L^Vi/i2K)j 



(B.J 



(B.9) 



4AVl^/(2/i) 

Note that both formulae reach the limit of pure diffusive 
transport for vanishing (compared to the diffusion coeffi- 
cient) values of the wind : 



\ pure — diffusion 



2K 



(B.IO) 



Appendix C. p and d fluxes from sources in the 
diffusive halo 

We now use a constant source term in the Galaxy: 

" ' (Vgai^)' + nd{z).v.a xN = q. (C.l) 



The general solution of this equation is given in Gradshteyn and Ryzhik"-^-^ 
(1994) (p. 1067, Eq. 9.255). Rewritten using z gives 

N{z) = Ae^"''^^ + Be^''''>\Tf{l^z). ^ function erf is defined as eri{z) 



2 

7^ 



dt. 



C.l. Constant wind Vgia{z) — Vc 

We proceed as for the standard source case. In the halo, 
the transport equation reads 

-KN" + VcN' = q . 

Using the boundary condition N{z ~ L) = 0, gives 

Integration on the thin-disk gives the relation 

-2KN'{0) + 2VcN{0) + 2hnv(jN{Q) = , 
so that the final solution is 



(C.2) 



iV(0) 



where 



qL ( 1 — (1 + a)e 



Ua + ^)(l-e-") + ae-" 



VcL hnvaL 
a^— and C = 



(C.3) 



C.2. Linear wind Vga.\{z) — Vi x z 

In the halo, the equation is 

-KN" + VizN' + ViN = q 



Using V = ^JVi/{2K) and following closely the derivation 
of Sec. B.1.2, we find the solution 

N{z) = Ae'"^^'' + Se(''")'erf(i.z) + -J. 
The condition N{z = L) = yields 

A=-Bx eTi[vL) - ^e-^''^'^\ 
Vi 

Integration of Eq. (C.l) over the thin-disk yields 

-2A'iV'(0) -f- 2hnv(jN{Q) = , 
which determines A and gives the final expression 



^ ' Vi\ e-t"^^) 



1 



where 



' — — and f 
2K ^ 



erf (i/z) — erf (I'i) 

^^/tt erf (z^L) 
2iyL + erf (i/L) 

hnvaL 



(C.4) 
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Fig. D.l. Exotic p flux for Kq = 1.35 X 10^^ cm^ s'^ (and 5 = 0.7), 
L = 4 kpc and Vc = 12 km s~^. Symbols correspond to the ID 
model with (plus) or without (open plus) spallative destruction. This 
is shown for an observer in the disk (2 = 0) or at ^ = 2 kpc. 

C.3. Formulae for = (no spallations) 



Constant wind Vr 



iV(0) 



qK 
V} 



1 




Linear wind Vi 



■1 iKy^ z 



iV(0) = -1(1 



(C.7) 



Finally, at sufficiently high energy, the purely diffusive 
regime is reached, leading, for both wind configurations, to 



N{z) 

N{Q) = 
^ ^ 2K 



q{L^ - z^) 



2K 



(C.8) 



Appendix D. The impact of neglecting spallations 

Solutions switching-on/off spallations for the constant 
wind model are compared using the master Eq. (C.2). 
Figure D.l displays a comparison of p fluxes for the best fit 
propagation parameters. The inclusion of spallations for p 
is of little consequence. We checked that it remains the case 
for other propagation parameters compatible with B/C 
data. The most important effect is when Vc = (which 
is usually never met). In that case, the maximum devia- 
tion is observed at low energy and is about several tens of 
percent — but we remind that due to the Solar modulation, 
fluxes below ^ 500 MeV/n IS energies are not reachable. 
The deviation is expected to be slightly larger for d, whose 
destruction cross section is roughly twice the p one. For 
consistency, the ID results were also compared with 2D 

calculations (sec companion) using qDarkiTi z) = \. 
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